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ABSTRACT 

An extremely sensitive all-sky survey will be carried out in the millimc- 
tre/submillimetre waveband by the forthcoming ESA mission Planck Surveyor. The 
main scientific goal of the mission is to make very accurate measurements of the spa- 
tial power spectrum of primordial anisotropies in the cosmic microwave background 
radiation; however, hundreds of thousands of distant dusty galaxies and quasars will 
also be detected. These sources are much more likely to be gravitationally lensed by in- 
tervening galaxies as compared with sources discovered in surveys in other wavebands. 
Here the number of lenses expected in the survey is estimated, and techniques for dis- 
criminating between lensed and unlensed sources are discussed. A practical strategy 
for this discrimination is presented, based on exploiting the remarkable sensitivity and 
resolving power of large ground-based millimetre/submillimctre-wave interferometer 
arrays. More than a thousand gravitational lenses could be detected: a sample that 
would be an extremely valuable resource in observational cosmology. 

Key words: methods: observational - galaxies: evolution - cosmic microwave back- 
ground - cosmology: observations - gravitational lensing - radio continuum: galaxies 



1 INTRODUCTION 

The properties of the population and the individual appear- 
ance of distant galaxies can both be modified significantly by 
the gravitational lensing effect of foreground masses. Grav- 
itational lensing has been investigated in most detail in the 
near- infrared/optical and radio wavebands, in which faint 
distant galaxies can be observed with sub-arcsecond resolu- 
tion. However, there are excellent prospects for extending 
these studies into the submillimetre waveband (Blain 1996, 
1997a,b, 1998). The fraction of gravitational lenses expected 
in a sample of galaxies selected in this waveband can be up 
to three orders of magnitude larger as compared with a sam- 
ple selected in other wavebands. 

The concept of magnification bias provides a useful way 
of describing the effects of lensing on a population of distant 
galaxies (Borgeest, Linde & Refsdal 1991). If the surface 
density of galaxies with flux densities greater than S v per 
unit redshift is n(Sv, z) at redshift z, then magnification by 
a factor A due to gravitational lensing would predict a mod- 
ified count, n'(S u ,z) = A" 1 nlSvA" 1 , z]. If n oc S" locally, 
then the magnification bias n'/n is given by A- (1+Q) . Hence, 
if a < — 1 then the surface density of galaxies is increased 
and the magnification bias is positive; if a > — 1 then the 
surface density of galaxies is decreased and the magnifica- 
tion bias is negative. Magnification bias was first discussed 
in the context of samples of bright quasars; however, counts 
of galaxies at faint flux densities in the submillimetre wave- 
band are expected to be uniquely steep, with a — —3 or 



less (Blain & Longair 1993, 1996), and so the magnification 
bias for faint galaxies in this waveband is expected to be 
very significant (Blain 1996). At brighter flux densities the 
submillimetre- wave counts are expected to be less steep, and 
to follow the Euclidean slope, with a — —3/2. The magnifi- 
cation bias is expected to be largest at flux densities at which 
the counts begin to rise above the Euclidean slope. If the 
population of dusty star-forming galaxies evolves strongly, 
as the first detections of previously unknown galaxies in the 
submillimetre waveband appear to demonstrate (Smail, Ivi- 
son & Blain 1997 - SIB), then the largest magnification 
bias is predicted to occur at a flux density of 0.1 to Uy, 
which is comparable with the sensitivity of Planck Surveyor 
(Bersanelli et al. 1996). 

The submillimetre-wave counts of lensed galaxies are 
expected to be sensitive to changes in both world models 
(Blain 1998) and scenarios of galaxy formation and evolu- 
tion. Hence, the properties of the lensed galaxies or quasars 
detected in the Planck Surveyor mission would allow the 
values of cosmological parameters to be investigated, even 
without detailed lens modeling. An Einstein-de Sitter world 
model and Hubble's constant Ho = SOkms" 1 Mpc" 1 are as- 
sumed throughout. 



2 THE PROPERTIES OF LENSES 

The surface and flux densities of gravitational lenses in the 
millimetre/submillimetre wavebands can be predicted by 
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Figure 1. The probabilities of galaxy— galaxy lensing as a func- 
tion of redshift assumed in this paper. The lensing optical depth 
to lensing with a magnification greater than two out to reshift z 
is a(z)/8. For references see Blain (1996). 



combining models of the population of distant dusty star- 
forming galaxies (Blain & Longair 1996) with a model of 
the magnification distribution due to lensing as a function of 
redshift (Peacock 1982; Pei 1993). The magnification distri- 
bution can be derived from the mass distribution of galaxies, 
and takes the form a{z)A~ i if A is large. Estimates of a(z) 
are presented in Fig. 1 for both evolving and non-evolving 
models of the distribution of lensing masses. In the evolv- 
ing model the mass distribution of the population of lens- 
ing galaxies is derived using the Press-Schechter formalism 
for structure formation by hierarchical clustering (Press & 
Schechter 1974), in which galaxies typically become smaller 
and more numerous at larger redshifts. The probability of 
lensing is predicted to be smaller in the evolving model as 
compared with the non-evolving model. The form of a(z) as- 
sumed in this letter is normalised to match the predictions 
of Pei (1993). 

Investigations of the evolution of galaxies and active 
galactic nuclei (AGN) in many different wavebands (Hewett, 
Foltz & Chaffee 1993; Dunlop & Peacock 1990; Oliver, 
Rowan-Robinson & Saunders 1992; Lilly et al. 1996) indicate 
that the observed evolution of both the global star-formation 
rate and the luminosity density of AGN is consistent with 
pure luminosity evolution (PLE) of the form (1 + z) 3 out to 
z ~ 2. The observations of SIB can be fit well by a model 
of PLE of the 60-/im luminosity function of IRAS galaxies 
(Saunders et al. 1990) with the form of (l + z) 3 out to z = 2.6, 
and a constant luminosity evolution factor of 46.7 in the in- 
terval 2.6 < z < 7. The predictions of this model, which is 
consistent with other recent observations by Kawara et al. 
(1997) and Wilner & Wright (1997) are shown in Fig. 2, for 
both an evolving and a non-evolving distribution of lenses. 

The effects of different world models are discussed by 
Blain (1997a, 1998) discusses the effects of different world 
models; a non-zero cosmological constant increases the pre- 
dicted surface density of lensed galaxies as compared with 
the Einstein-de Sitter model, because both the probability 
of lensing at a particular redshift is increased and a stronger 
form of evolution is required in order to explain SIB's ob- 
servations, as the smaller volume element at large redshifts. 
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Figure 2. The counts of lensed and unlensed field galaxies ex- 
pected with and without evolution of the lens population. 

The wavelength dependence of both the lensed and un- 
lensed counts is presented in Fig. 3. The principal feature 
is the large increase in the ratio of the surface densities of 
lensed and unlensed galaxies at flux densities of between 0.1 
and 1 Jy, which correspond to the onset of the steep rise 
in the unlensed counts above the Euclidean slope at these 
wavelengths. This increase is by almost two orders of mag- 
nitude as compared with the 175-/im counts, which have a 
form most similar to counts expected in the optical or radio 
wavebands. Another notable feature is the predicted increase 
in the surface density of unlensed 0.5-Jy galaxies as the ob- 
serving wavelength decreases. This increase is predicted to 
be larger by a factor of about 15 between wavelengths of 
850 and 350 /jm. The counts of lensed 0.5-Jy galaxies also 
increase with decreasing wavelength, from 1 mm to about 
200 fim; however, at wavelengths shorter than about 200 fim 
this effect is reversed, and the counts are predicted to de- 
crease as the observing wavelength decreases. 



3 A LENS SURVEY WITH PLANCK SURVEYOR 

Estimates of the numbers of galaxies that could be detected 
in an all-sky survey are listed in Table 1, based on the sen- 
sitivities quoted by Bersanelli et al. (1996) and the counts 
discussed above. In order to produce conservative estimates, 
an evolving lens population is assumed. Note that the pre- 
dicted number of lenses would be increased by a factor of 
about 5 if both a non-evolving lens model and a non-zero 
cosmological constant were assumed (Blain 1998). 

Between about 0.6 and 5% of the point sources detected 
by Planck Surveyor are expected to be lensed by a fore- 
ground galaxy. Despite a relatively large ratio of lensed to 
unlensed galaxies at wavelengths of 2000, 1380 and 850 /im, 
the surface density of detectable sources is expected to be 
relatively small at these wavelengths, as compared with the 
shorter wavelength passbands at 350 and 550 fim. The most 
useful wavelengths for detecting a sample of lenses are hence 
350 and 550 /im. More than a hundred lenses could be de- 
tected per unit solid angle in this systematic survey, and so 
even if about a third of these lenses are masked by Galac- 
tic emission at Galactic latitudes less than about 20 degrees, 
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Table 1. The sensitivities, source confusion limits and number of detections N expected in a Planck Surveyor survey. <7 c ; rrus and cTs are 
the expected Galactic cirrus confusion noise (Helou & Beichman 1990) and survey sensitivity (Bersanelli et al. 1996) respectively. cr c i rrus 
is calculated assuming that the mean Galactic background intensity at a wavelength of 100 (im is 5MJysr . S con f is the flux density 
at which the surface density of galaxies is expected to exceed 0.03 beam -1 ; the expected point source confusion noise <r con f ~ S con {/3, 
see Blain, Ivison & Smail (1998). An evolving population of lensing galaxies is assumed. In a map with Gaussian noise, about 10 5 , 2 
and pixels in the survey at 350, 550 and 850/^m would be expected to exceed flux densities of 3, 5 and 10<r due to random fluctuations. 
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many hundreds of lenses could still be detected in the sur- 
vey. This would increase the number of known galaxy-galaxy 
lenses by an order of magnitude. The bright counts of un- 
lensed galaxies in the submillimetre/far-infrared waveband 
would also be determined very accurately in such a survey, 
and so tight limits could be imposed to the form of evolution 
of the global star-formation rate at moderate redshifts. 



4 THE DISCRIMINATION OF LENSED AND 
UNLENSED GALAXIES 

How efficiently could a thousand lensed sources be sep- 
arated from a hundred thousand unlensed galaxies de- 
tected at 5a significance in the survey? Planck Surveyor 
will execute a multi-band survey, and so the colours of 
the detected sources could probably be used to reduce 
the size of the sample by a small factor. However, none 
of the sources will be resolved, and their positions will 
only be known to an accuracy of about ±2 arcmin. Other 
planned facilities - ESA's Far-Infrared and Submillimetre 
Space Telescope (FIRST) (Pilbratt 1997) and large ground- 
based millimetre/submillimetre-wave interferometer arrays 
(MIAs; Brown 1986; Downes 1996) - will allow this task to 
be completed in a practical amount of observing time. A 
strategy for following up the Planck Surveyor point source 
catalogue is outlined below. 



4.1 Pre-selection using Planck Surveyor colours 

Colour-colour and colour-magnitude diagrams for the 
sources detected in a simulated 100-deg 2 sub-field of the sur- 
vey are presented in Fig. 4. The distribution of lensed and 
unlensed galaxies are clearly different in each plot: lensed 
galaxies are found at large redshifts, and if the dust temper- 
ature in star-forming galaxies is correlated with their lumi- 
nosity, then lenses would be expected to have redder colours 
as compared with unlensed galaxies at the same flux densi- 
ties. The field of points in the colour-colour plot (Fig. 4a) is 
bounded at small and large redshifts by lines with a gradi- 
ent of about 0.6, which reflects the relative wavelengths of 
the three relevant observing bands. If the emissivity of dust 
grains was independent of wavelength, then both lensed and 
unlensed galaxies would lie on a single line with this slope, 
at positions determined by their redshift and dust tempera- 
ture. The points in Fig. 4(a) are spread within a box-like re- 
gion because a wavelength-dependent emissivity is assumed. 
These plots can be used to determine which sources are 
probably at low redshifts and to prioritize the sample to en- 
hance the probability of detecting lensed sources earlier in 
the follow-up work; the subsequent follow-up stages should 
proceed through the sample working from the bright red 
sources to the faint blue ones. 



© 0000 RAS, MNRAS 000, 000-000 



4 A. W. Blain 

(a) (b) 




1 2 0.1 0.2 

S 350 / S 550 S 550 / J ^ 



Figure 4. Colour— colour and colour— magnitude diagrams for both lenscd and unlensed galaxies in a simulated 100-deg 2 sub-field 
of the Planck Surveyor survey. The flux densities of the plotted galaxies exceed the 5cr 550-/iin sensitivity of the survey; in (b) they 
also exceed the 5(7 sensitivity at 350-/im. Lensed galaxies are always found at rcdshifts greater than unity, and are typically redder 
as compared with unlensed galaxies. 



4.2 Determining more accurate positions 

FIRST will be equipped with a submillimetre-wave bolome- 
ter array receiver with a 6-arcmin field of view and a res- 
olution of about 20arcsec (Griffin 1997). This field of view 
is accurately matched to the resolution of the pixels in the 
Planck Surveyor map. The 350-/im flux density of candi- 
date lenses is expected to exceed about 0.2 Jy, Fig. 4(b), 
and so each source could be detected at 5a significance and 
located by FIRST in a 20-s integration. Assuming 200 per 
cent overheads, follow-up observations of 10 5 lens candidates 
would take about 70 days of integration, a period compa- 
rable with the time required for a shallow FIRST galaxy 
survey (Rowan- Robinson 1997). A ground-based instrument 
such as SCUBA at the JCMT (Cunningham et al. 1994) 
would require four separate pointings and a total integra- 
tion time of about 3 hours at 850 /im in order to detect a 
candidate lens at 5a significance. Hence existing ground- 
based facilities are not sufficiently sensitive to carry out this 
programme. 

A large MIA is expected to detect a 50-mJy 850-/im 
source at 5a significance in about 28 ms, but can only image 
about 0.08 arcmin 2 of sky per pointing. Hence about 250 
separate pointings, and a total integration time of about 7 s 
would be required to image the field of each Planck Surveyor 
pixel that contains a candidate lens with such an instrument. 

Either a space-borne submillimetre-wave telescope or 
a ground-based submillimetre-wave interferometer could lo- 
cate all the potential lensed sources in the Planck Surveyor 
catalogue to an accuracy better than about 20arcsec in 
about 3 months of observing time. Observations with FIRST 
have the advantage of simultaneously determining the spec- 
tral energy distributions of the candidates at wavelengths 
shorter than 450 /xm. 



4.3 Diagnosis of lenses 

Once reasonably accurate positions are known for the candi- 
dates, an MIA can be used to image them at a large signal to 
noise ratio and search for signs of lensed arcs, rings or mul- 
tiple images. The relative numbers of 0.1- Jy lenses per unit 



magnification are expected to be flat out to magnifications 
of about 50 (Blain 1997a; Fig. 4), and so most detected lens 
should show these clear signatures of strong lensing. The 
flux densities and redshifts of the candidates are expected 
to be similar to that of the z = 2.8 starburst/ AGN recently 
discovered by Ivison et al. (1998). Several hours of integra- 
tion were required to image and take the spectrum of this 
source using 4-m class optical telescopes, and so even with 8- 
m telescopes available it will be quite impossible to follow-up 
the Planck Surveyor catalogue in the optical/near-infrared 
waveband in a reasonable time. 

In a 1-minute 850-fim integration on a 2-km baseline, 
an MIA could reach a 5a flux density limit of about 1 mjy 
at sub-O.l-arcsecond angular resolution. This corresponds 
to the detection of a lens candidate in the Planck Surveyor 
catalogue at a significance of about 250. Any multiple im- 
ages, arcs or rings in the source will be clearly resolved us- 
ing such an instrument. Imaging 10 sources in this man- 
ner would take about 70 days of dedicated observations. Of 
course many of the candidates will show no signs of multiple 
structure, and so could be excluded from the list of candi- 
dates in a shorter integration. Note that the lensing galaxy 
is expected to be optically thin to submillimetre-wave radia- 
tion and so differential extinction of multiple images, which 
affects lens surveys in the optical waveband, will not compli- 
cate the selection effects in a Planck Surveyor survey. Pos- 
sible other selection effects are considered in the following 
section. 

Ultimately the catalogue of good candidates derived 
from the high-resolution MIA images will require spectro- 
scopic confirmation in the optical/near- infrared in order to 
diagnose them as lenses unequivocally. The results of the 
large-area Sloan Digital Sky Survey (Kent 1994) will prob- 
ably be of great assistance in identifying potential lensing 
galaxies. 



4.4 Potential problems and selection effects 

Chance superpositions of sources, multiple components of 
an unlensed source and emission from the lensing galaxy 
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could bias the selection of candidates in the Planck Surveyor 
survey and confuse the diagnosis of lenses. 

First, the counts at flux densities comparable to the 0.1- 
Jy selection limit are expected to be about 3deg -2 , Fig. 2, 
and so the chance alignment of two sources of comparable 
brightness within about 5arcsec of each other, which could 
be mistaken for different components of a multiple image, 
is expected to be of the order 0.001 per cent, equivalent to 
less than one source in the whole catalogue. About 1.5 per 
cent of the pixels in the Planck Surveyor survey are expected 
to contain more than one detectable source; however, these 
sources will typically be well separated in follow-up obser- 
vations at finer angular resolution. 

Secondly, although interacting galaxies show clearly 
separate components in the optical waveband, high- 
resolution millimetre- wave observations show that dust radi- 
ation typically originates from a single core of order several 
hundred parsecs in size (Solomon et al. 1997). At z = 2 this 
corresponds to an angular scale of a few times 10 -2 arcsec, 
a point source for the purposes of this investigation. Hence, 
because Planck Surveyor will detect only the most luminous 
distant galaxies, it is unlikely that intrinsic multiple struc- 
ture in distant galaxies will lead to the false diagnosis of a 
large number of lenses. 

Thirdly, the emission from a lensing galaxy would make 
the detection of a lens candidate more likely in the Planck 
Surveyor survey, and complicate its diagnosis. However, al- 
though in the optical waveband lensed images are typically 
fainter than the lensing galaxy, the situation in the submil- 
limetre waveband is much more similar to that in the radio 
waveband, in which the lensed images are much brighter 
than the lensing galaxy. This effect was discussed by Blain 
(1997b) in the context of confusion between distant lensed 
submillimetre-wave galaxies in the field of a cluster and 
the cluster galaxies themselves. In Fig. 5 the flux density- 
redshift relations expected for an V IRAS galaxy, with a 
luminosity of about 6 x 10 L0 at the present epoch, and a 
detectable lens candidate, assumed to be at z > 0.5 are com- 
pared. Unless the lensing galaxy is at z < 0.1, which seems 
rather unlikely for the typical lens candidate at z » 1, the 
flux density from the lensed images is expected to domi- 
nate that of the lens, even if the lensing galaxy was a lu- 
minous IRAS galaxy. Most lenses would be expected to be 
quiescent elliptical galaxies with much smaller far-infrared 
luminosities. If a severe colour-based pre-selection was im- 
posed, then perhaps contamination by the lensing galaxy 
could introduce a small bias into the sample; however, we 
have shown that the follow-up of the Planck Surveyor survey 
is still practical even without any pre-selection. 



5 THE OPTIMAL STRATEGY 

It will be practical to compile a sample of about a thou- 
sand lenses, starting from the 10 5a detections expected in 
the Planck Surveyor point source catalogue, by conducting 
a concerted programme of observations using FIRST, large 
MIAs and conventional optical telescopes lasting several 
months. The most promising candidates can be pre-selected 
for follow-up observations based on their submillimetre-wave 
colours in the catalogue. FIRST or an MIA can then be 
used to determine the positions of the candidates with suf- 
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Figure 5. The flux density— redshift relations expected for an 
L* IRAS galaxy and a lensed object detected at her at 550 /im 
in a Planck Surveyor survey. A low dust temperature of 30 K is 
assumed to boost the low-redshift flux densities. The upper and 
lower dotted lines correspond to the approximate lcr sensitivities 
of the Planck Survey or survey, Table 1, and a 1-minute integration 
with a large MIA respectively. 



ficient accuracy to permit a programme of high-resolution 
follow-up imaging observations using an MIA. Any lensed 
arcs, rings and multiple images in the candidate sources will 
be identified at this stage. Finally, the remaining sample of 
about a thousand excellent candidates would be subjected 
to spectroscopic confirmation or refutation as lenses using 
8-m class telescopes in the optical/near-infrared wavebands. 



6 CONCLUSIONS 

The Planck Surveyor mission has the potential to detect 
a very large sample of gravitationally lensed galaxies and 
quasars. The exact size of the sample will depend on the 
properties of distant dusty galaxies and quasars, the cosmo- 
logical parameters, the observing wavelength and the sensi- 
tivity of the survey, but about a thousand lenses and about 
a hundred thousand unlensed galaxies are expected to be de- 
tected a 5a significance at wavelengths of 350 and 550 /im. A 
practical strategy for submillimetre-based follow-up obser- 
vations has been outlined using FIRST and a large ground- 
based millimetre/submillimetre-wave interferometer array. 
The scientific rewards from the compilation of such a large 
unbiased sample of lenses would be very considerable, and 
it is difficult to see how such a sample could be compiled in 
any other waveband. 
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